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Abstract 

In a continuation of an ongoing program, we use staggered chiral perturbation theory to calcu- 
late the one-loop chiral logarithms and analytic terms in the pseudoscalar meson leptonic decay 
constants, f w + and f K +- We consider the partially quenched, "full QCD" (with three dynamical 
flavors), and quenched cases. 
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I. INTRODUCTION 



Simulations with staggered (Kogut-Susskind, KS) fermions are very fast relative to other 
available approaches, making possible simulations of QCD that include the effects of light 
sea quarks ll | . However, with currently practical lattice spacings (e.g., MILC simulations 
2I S, 0, S, l| at a ~ 0.09 — 0.13 fm) taste 1 violations are not negligible. Thus fits to such 
lattice data should take into account the taste-violating effects; indeed, if such effects are 
not taken into account, the speed advantage of KS fermions may be offset by the size of the 
systematic errors. The taste-violating effects can be calculated in a systematic way using 
staggered chiral perturbation theory (SXPT). 

In Ref. we formulate SXPT for the physical case of multiple flavors. SXPT is then 
used to calculate the one- loop chiral logarithms in the pion and kaon masses. Here, we 
continue the program of Ref. (7[ and compute f n + and f K + , the tt + and K + leptonic decay 
constants for the Goldstone mesons, to one loop. As we have laid most of the necessary 
groundwork already, we will merely state what is necessary for this present work and refer 
the reader to Ref. pj for the details common to both calculations. As in the calculation of 
the 7r + and K + masses, we perform our calculation using three dynamical KS-flavors (each 
with four tastes), which we call the 4+4+4 theory, and later adjust the result by hand using 
a quark flow technique 2IEI to a 1 + 1 + 1 theory (three flavors each with a single taste). 

The outline of this paper is as follows: In Sec. |Hj we write down the SXPT Lagrangian 
for three dynamical flavors. We then calculate, in Sec. 1111} the one-loop chiral logarithms 
which contribute to the flavor-nonsinglet Goldstone meson decay constant in the partially 
quenched case. Here we keep three dynamical flavors but add two additional quenched 
flavors as valence quarks, which in the general case have distinct masses from the dynamical 
(sea) quarks. The transition to a 1+1+1 theory is then made. There are only a few differences 



in this procedure from that of Ref. [7]. We also give the results in the fully quenched case. 
The full next-to-leading order (NLO) results, including the analytic terms, are presented 
in Sec. IIVI for various relevant cases. We conclude with some comments in Sec. An 
Appendix gives technical details about the evaluation of the one-loop integrals that arise in 

Sec. urn 



1 We use P, y| the term "taste" to denote the different KS species resulting from doubling, and "flavor" for 
the physical u-d-s quantum number. 
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II. THE LEE-SHARPE LAGRANGIAN FOR 3 FLAVORS 

The starting point for SXPT is the Lee-Sharpe Lagrangian |icj generalized to multiple 
flavors. In Ref. Q we examined a general n-flavor theory* and later specialized to „ = 3. 
Here we take n = 3 from the beginning. For 3 KS flavors, S = exp(z$//) is a 12 x 12 
matrix, with $ given by: 

f U 7T+ K+\ 

$ = 



V 



u 


7T + 


K+ 




D 


K° 


K~ 


K° 


S 



(1) 



/ 



where U = Yl a =i Ua,T a (and similarly for ir + , K + , . . . ), with 

T a = {£5, i£p5, i^u, £m £/}■ 



(2) 



We use the Euclidean gamma matrices with ^ u = (// < v in Eq. (J2J)), ^5 = £^£5, 
and £/ = / is the 4x4 identity matrix. The field S transforms under S{7(12)i x S , C/(12)^ 
as S — > LHw . The component fields of the diagonal (flavor-neutral) elements (?7 a , _D a , and 
S a ) are real; the other, charged, fields (7r+, etc.) are complex, so that $ is Hermitian. 
The mass matrix is given by the 12 x 12 matrix 

/ ™ t n n \ 



M 



mj 
m d I 
mj 



(3) 



Our (Euclidean) Lagrangian is: 



f 2 1 9tt? 2 

^-Tr(^E«9 M £t) _ ^TrCME + TW^t) + ^ Uj + Dl + gtf + ^ (4) 
o 4 o 



where /1 is a constant with units of mass, Tr is the full 12 x 12 trace, and V = 14 + 14' is 
the taste-s—v breaking potent... The M ter,„ is g,ve„ in Ref. Q it is not needed 



Here n refers to the number of sea quarks. 
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explicitly here. For U we have 

-W = C 2y ^[Tr(^ 3) S)Tr(^ 3) S) + / i .c.] 

V 

+C 2 A^[Tr(^ 3 5 ) S)Tr(^ ) S) + / i .c.] 

V 
V 

+C5A^[Tr(ei 3 5 ) S)Tr(e£ ) St)], (5) 



2 

V 



(3) 

where the £ B are block- diagonal 12 x 12 matrices: 



t(3) 
SB 



(6) 



(& o o \ 

o 6, o , 

V (J 

with £ B the 4x4 objects, and B e {5, //, /iz/ (/i < z/), /x5, /}. 

As seen in Ref. jjj, U' generates two-point vertices at 0(a 2 ) (shown in Fig. QJ that mix 
flavor-neutral particles of vector and axial tastes. In addition, flavor-neutral, taste-singlet 
particles are mixed by the wig term in Eq. (J1J, which results from the anomaly. In all three 
cases (taste vector, axial vector, and singlet), we have a term in the Lagrangian of the form 
(5'/2)(U + D + S) 2 , where 

{a 2 S' v , taste- vector; 
a 2 5' A , taste-axial; (7) 
Aml/3, taste-singlet. 

Expressions for S' v and 8' A in terms of the coefficients of U ' are given in Ref. |7j . These 
mixings require us to diagonalize the full mass matrix in each of the three channels. We 
write the propagator for the vectors as: 

G v = Goy + V v . (8) 

D v is the part of the taste- vector flavor-neutral propagator that is disconnected at the quark 
level (i.e., Fig. [T]plus iterations of intermediate sea quark loops). Explicitly, we have [2J: 

v v _ 2« (g 2 + m 2 Uv ){q 2 + m 2 Dv ){q 2 + m 2 Sv ) 

MN ~ a V {Q 2 + rn 2 Mv ) (<7 2 + < v ) (<7 2 + ™ 2 o ) (d 2 + < ) (<1 2 + ™\ ) ' K ) 
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Here, m^ , m 2 v and m 2 ^ are the eigenvalues of the full mass-squared matrix (i.e., the poles 
of GV). We emphasize that Eq. @ remains valid in the n = 3 partially quenched case. The 
external mesons M and N may be any flavor-neutral states, made from either sea quarks or 
valence quarks. 

In the quenched case T>)^ N is simply 

MN ~ V tf + m* Mv W + m% v ) ■ (10) 
Equations (jHJ) through (JTUjl apply explicitly to the taste-vector channel; to get the result 
in the taste-axial (taste-singlet) channel, just let V — > A (V — > I and a 2 5' v — > 4mg/3). 
In the quenched case we cannot take mo — > oo, and must include additional ?7j-dependent 
terms in the Lagrangian, resulting in the replacement — ► ml + aq 2 in the singlet form of 

Eq. jnu 0,0. 



III. ONE LOOP DECAY CONSTANT FOR 4+4+4 DYNAMICAL FLAVORS 



We calculate the pion 3 decay constant in a partially quenched theory. Full theory results 
are easily obtained by taking appropriate limits. There are three sea quarks (u, d and s) 
and two valence quarks (x and y). The pion of interest is the P 5 + , a Goldstone pion which 
is composed of an xy pair of quarks. 

The P 5 + decay constant is defined by the matrix element: 







J fi5 



P 5 + (p)) = -z/p + p M 



(11) 



where is the axial current corresponding to P 5 + . With this normalization, f n w 131 MeV 
In terms of E, we can write the axial current as 



J fib 



8 



-Tr 



(12) 



Here V p+ projects out the 4x4 block with appropriate flavor: If we make x and y the last 
two flavors of S, then = ^5^4, where i, j are flavor indices. 
At one loop, the decay constant has the form 

1 



/r+ = / 1 + 



16tt 2 / 



2 f2 J P; 



(13) 



We refer generically to any flavor-charged meson as a "pion." 
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There are two contributions to Sf P +, which we call Sf^ + and 5f^" ent . They are shown in 
Figs. El and El respectively. The contribution 5f^ + is merely wavefunction renormalization. 
We have: 

z I 167T 2 / 2 d^) 

Sf P + = = 2- ■ (14) 

The self-energy, X(p 2 ), has already been calculated in Ref. pj. The wavefunction renormal- 
ization arises only from the vertex generated by the kinetic energy term in Eq. (J3J), since 
derivatives on the external lines are necessary to generate p 2 dependence in a tadpole dia- 
gram. The factor of 1/2 in Eq. (|14jl is due to the fact that this diagram is multiplied by 



*Z = Vl + SZ « 1 + \5Z. 
The contribution 8f™* Tent is the current correction. It arises from the expansion of Eq. (|12|) 

5 

to cubic order in Performing this expansion, it is easy to see that the wavefunction and 
current correction terms are proportional to each other: 8f^+" axt — —4Sfp+. This fact, 
noted also in [9|, is perhaps not surprising, since the form of the axial current, Eq. (|12p. is 
determined through Noether's theorem only by the kinetic energy part of the Lagrangian. 
From Eq. ()14j) . we then have 

5f P + = %l + + 5fpT ent = -\SZ p+ . (15) 
Using intermediate expressions for X(p 2 ) from 7], the one- loop result is 



% = -i / 5 



7T 



1 



Q,B H 



(16) 



Here, Q runs over the six mesons formed from one valence and one sea quark {i.e., the xu, xd, 
xs, yu, yd, and ys mesons). As before, B takes on the 16 values {5, /i, \iv{p, < v),jih, I}. We 
have already included the factor of 4 that comes from summing over the degenerate vector 
and axial contributions in the D v and T> A terms. Despite the fact that the only 4-point 
vertices contributing to this expression come from the kinetic energy term, the result is more 
complicated than that for the mass renormalization Q because there are no cancellations 
here (either accidental or required by symmetry). 

The first term Eq. (|16J) (the sum over Q and B) comes from the wave function renormal- 
ization and current correction diagrams shown in Fig. Ufa), which involve a single virtual 
quark loop. The diagrams arise from the vertices in Fig.|^a), respectively, where i is summed 



over the sea quarks only. In this case, the propagator in the loop must be connected since 
the loop meson is not flavor-neutral. 

The vertices in Fig. EJa) also produce diagrams with disconnected loop propagators, 
Fig. E(b) and (c), when i — y (or i = x in the y <-> x version of Fig. E[a)). These diagrams 
give rise to the V Y y and V X x terms in Eq. (jTfij) . 

Finally, the vertices in Fig.|3fb) generate the diagrams Fig. QJd) and (e). The T>xy terms 
in Eq. (fTBj) come from these diagrams. For more discussion of how to identify quark flow 
diagrams with the SXPT contributions, see Ref. 

We can write down the quenched result by (1) eliminating the term summed over Q 
and B, which arises from virtual quark loops (diagrams Fig. 11(a)), and (2) replacing of 
T> — > X>i uench for T> v , T> , and T* 1 . These replacements eliminate diagrams Fig. 0(c) and (e). 

In the partially quenched case, when the x or y quark mass is different from all sea 
quark masses, there are double poles in Eq. (fiTjj) coming from the T>xx or Dyy terms. This 
is different from the mass renormalization result, where double poles do not arise unless 
m x = m y and this mass is different from all sea quark masses, a case we did not treat in 
detail in Ref. [?J. In order to write down explicit results for partial quenching here, we must 
therefore expand on the notation of Ref. [7] . Below we will use the notation defined in the 
Appendix, where we explain the techniques we use to evaluate the integrals. 

Before performing the momentum integrals, we now make the transition from the 4+4+4 
case to the 1 + 1 + 1 case. This is easily accomplished since we have already determined the 
separate contributions from diagrams in Fig. H] with various numbers of sea quark loops. 
Those diagrams with a connected propagator in the loop, Fig. 11(a), have a single sea quark 
loop and simply must be divided by 4. 

The remaining diagrams have the same form as those treated in Ref. |t], so we just 
briefly review the procedure. Diagrams (b) and (d) in Fig. H] have no sea quark loops and 
are unaffected by the transition to the 1 + 1 + 1 case. These diagrams have a single factor of 
5' coming from the overall coefficient of the disconnected propagator T> in Eq. (JOJ) . Each sea 
quark loop added on to diagrams (b) and (d), as in (c) and (e), comes with an additional 
factor of 5'. Therefore, we must merely make the replacement 5' — > 5' /4 in all but the overall 
factors of 8' . This is easily accomplished by letting 5' — > <5'/4 in the computation of the full 
mass eigenstates (i.e., m^ , m^ v and m^, in Eq. ©), but not in the overall coefficient of T>. 

After making the transition to the 1+1+1 case, taking the m — > oo limit (with m^, ~ wIq), 



and using Eq. (|3ip through Eq. (j4"Uj) to perform the momentum integrals, we have 

sf Ps+ - ~ KJ + ^fe 31 ({-MS})^) + 4f({M^})£(m 2 Yi ) 

Q,B ^ 

+ E *%2t,i{M$)mo + E ^((KwK) 



.7/ 



J I 



.7/ 



< 3] ({-MgjM™*v 



+ 



(17) 



+ £ ^S5v({^MO + 2 E 

where Q and S have the same meaning as in Eq. (JTBJ), the chiral logarithms £(m 2 ) are defined 
in Eqs. (|43|) and ()45|) for infinite and finite spatial volume, respectively, and the Rs and .Ds 
are residues defined in Eqs. (|33j) and (|36J1 . The arrow signifies that we are only keeping the 
chiral logarithm terms in this expression. We have defined the sets of masses in the residues: 

{Mz } = {m^o,m v ,m z } , 
{.M (2) } = {irv,ffl,,fflxiWy}i 

{.M (4) } = {171^,7X1^ m v >,mx,m Y } , (18) 

where Z can be either X or Y, and we show the taste labels explicitly in Eq. ()17j) . We do 
not include the numerator masses in the argument of the residues, as they are the same for 
each case: 

{//} = {mu,m D ,m s } , (19) 

with appropriate taste subscripts. The sums over jj, jy, and ja run over the set of masses 
included as the argument of the residues in each sum. 

The values of m 2 , m 2 and m 2 , in each taste channel in Eq. (fTTj) are the eigenvalues of 
the full mass matrix 



(m\ + 5' 74 S'/A 
5' I A m 2 D + 5' I A 



\ S'/A 



S'/A 



m\ 



S'/A \ 
5'/A 

+ S'/A/ 



(20) 



where S' is given by Eq. (J7|), the masses mfj, m 2 D , m| have an implicit taste label (V, A, or 
/) depending on the channel, and the m 2 , — > oo limit should be taken in the singlet channel. 



The explicit expressions for these eigenvalues are not illuminating in general, but they are 
given for m u = rrid (the 2 + 1 case) in Ref. Q. 



In writing down Eq. (|17|). we have assumed m 2 x ^ m Y . When m 



x 



m Y some of 



the residues of sets {A4^} and {.M 1 ' 4 )} blow up, so the limit must be taken carefully. 
Alternatively, one can simply return to Eq. take the limit trivially, and perform the 

integrations again. 

In the quenched case, we drop the first term in Eq. (j!6|) (with the sum over Q and B) 
since it comes from diagrams with a sea quark loop. In the remaining expression, we make 
the replacement V — > £> quench . 

Writing out the residues explicitly, and using the results from the Appendix, we obtain 



Sf 



quench 



771, 

~6 



[m 



Yj 



m Yi - m 2 Xi 



a 
+ 6 



1 

2 ( 



+ '-a 2 8' v 



i(m Xl )+£(m Yl ) 
i(m Xl ) - m Xi £(m Xi ) + £(m Yl ) - m Yl £(m Yl ) + 2 



in 



X, 



{m 2 Xi )-ml 



l Y v y 



[m 



ni 



Yi 



m 



Xy j 



rriv — m 2 Y 



+ 



V —> A 



(21) 



Carefully taking the limit m y 



m 7 



(or returning to Eq. ()16|) and taking the limit trivially), 



we see that the singlet terms vanish but the vector and axial terms do not. This is consistent 



with the known result 



ft 



12j in the symmetry (continuum) limit that there are no chiral 



logarithms in the quenched pion decay constant with degenerate masses. 



IV. FINAL NLO RESULTS 

For the complete expression for the NLO pion decay constant, we need the u O(p 4 )" 
analytic terms in addition to the chiral logarithms calculated above. Ours is a joint ex- 
pansion in a 2 and the quark mass m, so we are looking for the analytic contributions 
arising from terms of 0(m 2 ,ma 2 ,a 4: ) in the chiral Lagrangian. Examples of such terms 
are: Tr(d^d^)Tr(MZ + MTt) [0{m 2 )}, a 2 VTi{d^d^) [0(ma 2 )}, and (a 2 V) 2 [0(a% 
There will also be chiral representatives of those 0(a 2 ) operators in the effective continuum 
QCD action ("Symanzik action") that have no representatives at lowest order — such op- 
erators comprise what Lee and Sharpe call Sq F< " B \ Chiral representatives of Sq F ^ BS> 

I I 

operators have two derivatives 10] and therefore are 0(ma 2 ). 
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The only terms in the chiral Lagrangian that can contribute to the decay constant are 
those with derivatives, so 0(a 4 ) terms are not relevant here. Similarly, the "m" in 0(ma 2 ) 
terms must come from two derivatives. Therefore such terms just make a NLO contribution 
to f p + of the form Fa 2 f, where F is a constant formed out of the coefficients of the relevant 
Lagrangian terms. F would of course depend on the taste of the decaying particle, but we 
are considering only Goldstone particles here. Finally, the terms in the Lagrangian that are 
0(m 2 ) are just the NLO ones familiar from continuum XPT 



A. Full and partially quenched NLO results 

Using the definitions of Li in Ref. we thus get from Eqs. (fTTj) and ([T3]). in the 1+1+1 
partially quenched case, 

1 



./ 



l-loop, 1+1+1 



/I 



16vr 2 / 2 



Q,B V 

+R^ ] ({M^})£(ml l ) +J^D%%({M { i]})£(ml) 

ji 

+ E ^%({M^})£(ml) - 2 J2 R [ *?\{Mf})l{ml 



31 



+ £^:l, ( {>f«> })<(<) +£^U*0)«(- 

3v 



,(3) 



31 

2 x , D [4,3], 



,(3)- 



3V 

16/i 



3V 

+ (V ^ A 



P 



8u "1 

m u + m d + m s ) L 4 + — (m x + m y ) L 5 + a 2 Fj . 



(22) 



Definitions here are the same as in Eq. (|17|) . We have checked that this result reduces to 
that of Sharpe and Shoresh [ijj in the continuum (symmetry) limit. Using Eq. (jHEJ), it is not 
hard to show that changes here in the chiral scale A can be absorbed into the parameters 
L 4 , L 5 and F, as expected. 

In the 2+1 case (m u = ma = rug) with no other degeneracies, there is some simplification 
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because m^ = mfj = m? D in each taste channel. We obtain: 



^iccmw =/ l 1+ 1 



16tt 2 / 2 



2 



- 2 ^^({M 6) + ^ 2 ^fe 21 ({-M2U)^) 

h 

jv 

+T,^i^MO + »E^({^})<«)) + - -i; 

jV jv 

16/x 8/i 2 1 

+-T2" ^ 2m ^ + m ^ + J2 ( mx + -^5 + a F| , (23) 

with definitions the same as in Eq. (|17p. except that now the denominator masses in the 
residues are: 

{.M^} = {m v ,m z } , 
{.M (6) } = { mr} ,m x ,m Y } , 
{M?} = {m v ,m v ,,m z }, 

{M {8) } = {171^,711^, m x ,m Y } , (24) 

where Z can again be either X or Y, and a taste label is implicit. The numerator masses in 
the residues of Eq. (J23p are not shown explicitly. They are always 

M = {mu,m s } , (25) 

with the taste label again implicit. 

Various cases of interest can be obtained either by carefully taking limits in Eq. (j22j) or 
(I23j) . or by taking the limits in Eq. |T6|) and redoing the momentum integration. We first 
consider the "full QCD" case of "real" pions and kaons. By setting m x = m u and m y = ntd, 
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but keeping m u ^ rrid, we get after a bit of algebra: 

1 



l-loop, 1+1+1 



/I 



16tt 2 / 2 



{m% 1 - m 2 Di ) 2 



Q,B 

({M 2) };{™D>K) 



j i 



JV 



3V 



1 6 fi 'i i .-) 

+ —zr {m u + m d + m s ) L 4 + — (m u + m d ) L 5 + a F 



8/i 
T 2 



(26) 



where the sets {Ai^} and {.M'- 4 *'} are given in Eq. (fT%|) (with X — > U and Y — > D), jj and jv 
run over all masses in {M. }} and {A4 V }, respectively, and Q G {U, D, n + , n~ , K + , K }. Note 
that there are no double pole terms here, due to cancellations in the disconnected flavor- 
neutral propagator. The charged kaon result can be obtained from Eq. ()26|) by making 
the replacements d «-> s and D «-> S wherever they appear explicitly (as well as in the 
definitions of the mass sets {.Mf-} and {.My} and in the sum over Q, where we now have 
Q E {U,S,K+,K-,7c + ,K }). 

The full theory pion result simplifies even more in the 2 + 1 case, when the up and down 
quark masses are equal. After writing the residues explicitly, we obtain: 

1 

16 



l-loop,2+l 



/ 1 + 



16tt 2 / 



+2a 2 5 



2P 
V 



.m 



B 

2 

— Tfl 



b v v v 



in 



s v 



(m 2 Vv -ml 9 ){m 2 n , - ) 



'Iv 



i(mlo v )\ + (V 



A 



16/i 



16/i 



(2m^ + m s ) L 4 H — — m^L 5 + a F > . 



(27) 
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Similarly, the 2 + 1 result for the full kaon is: 



/ 



l-loop,2+l 



f 1 + 



16tt 2 / 2 



-l^(2£«o) + 3£K, 



+ £ mt 



+- 



K ? ) - -£(r<) + -£(m^ 



+ -a 2 5' v 



,m 



Sv 



+ m o — 2m 



Vv 



m 2 A - m^ v ){m 2 Sv - m^m 2 ^ - 

(m| + m 2 - 2m 2 , f 

_ ZX. V JL £( m ', ) 

(m 2 —m 2 , )(m% —m 2 , )(ml —m 2 , ) Vv 

K 7r?y rf v /\ S v rf v /\ Vv rfyl 

,2 



t{m. 



2 \ 

w 



w 

2 



in 



2 

Sv 



m\ 



<° ) + 



m\ 



m 



S v 



« - m \ ) ( m l' v ~ m l v ) ^ ( m vv - m lv ) - m lv ) 



-£(m 



2 ' 
Sv - 



+ V-> A 



16// 



8/x 



+ -— (2m £ + m s ) L 4 + — (m< + m s )L 5 + a F 
/ / 



(28) 



Here we have used the fact that m 2 ^ 
It is easy to check that Eqs. (J27j) and 
limit, where all tastes are degenerate. 



r + |^^o in the 2+1 case to simplify the result. 



reduce to the standard answers 



1311 in the a 2 —>■ 



B. Quenched NLO Results 



In the fully quenched case, we only need to consider the two cases m x ^ m y and m 3 
For the quenched "kaon" case (m x ^ m y ) we obtain: 



f 



1— loop, quench 
K+ 



/ I 



167T 2 / 2 



£(m 2 Yi ] 



in 



Yi 



m 



a 



+2 



m 2 x /(m 2 Xi ) - m 2 Yi £(m 2 Yi ) 



+ -a%[l(m Xv )+£(m Yv ) 



+2- 



(m 



(m 



Yv> 



m Y v ~ m X v 



+ V 



A 



+ j^(m x + my) L' 5 + a 2 F' \ . (29) 



The analytic terms in the quenched case are marked with primes to indicate that they 
may have different values than in the full theory. Also, note that there is no analytic term 
involving the sea quarks in the quenched they play no role here. In the continuum 

limit, Eq. (|2*9*|) reproduces the known quenched result [l2j. 
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Taking the degenerate limit (m y = m x ) in the quenched case, we obtain for the quenched 
"pion" : 

/ i-ioo P ,« h = / | 1 + _^_ ^a%I(m 2 Xv ) + 2a'5'J(m 2 XA j) + ^m x L' 5 + a 2 F'}. (30) 

This is consistent with the fact that in the isospin limit, the continuum quenched pion decay 
constant does not contain chiral logarithms. 



V. REMARKS AND CONCLUSIONS 



The most general result we have is for the n = 3 partially quenched case (1 + 1 + 1) with 
all valence and sea quark masses different, Eq. (J22j) . Other interesting cases can be obtained 
from Eq. ()22)) by taking appropriate mass limits. The results most relevant to current MILC 
simulations are those with m u = = mi (the 2+1 case); these and other important limits 
are presented explicitly in Sec. IIV Al The results in the quenched case are given separately 
in Sec. HVBl in Eqs. flU and 

The explicit results in Sec. II VI often appear dauntingly complex. However, the intricacies 
arise primarily from the momentum integration, which produces chiral logarithms with com- 
plicated residues from each of the many poles in the disconnected flavor neutral propagator, 
Eq. (jHJ). The result before integration, Eq. (fTtjj). is actually quite simple, and the reader may 
prefer to start with that expression and perform the integration himself in specific cases of 
interest . 

In the partially quenched case, double poles arise here even when the valence masses are 
non-degenerate, just as they do in the continuum [liiQl- It is interesting that these double 
poles appear in the explicitly 0(a 2 ) terms (taste-vector or axial channels, proportional to 
5y or 5' A ) as well as in the continuum-like taste-singlet channel. 

Using the SXPT results presented here and in [3], it seems possible to fit existing lattice 
data and extract physical physical parameters (e.g., f^, fx, fn s , {m u +md) /2, Li) with rather 



small discretization errors 



The next steps would be to extend the current approach to 



describe heavy-light mesons [12( and baryons. 
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APPENDIX 

In this Appendix, we go through the technical details of calculating the integrals found 
in Eq. For the terms only containing single poles, this was done in Ref. jQ], so here we 
focus on the terms which contain double poles. 

Consider first an integrand of the form 

2™ ({m};M) - , (31) 

where {m} and {//} are the sets of masses {mi,m 2 , . . . ,m n } and {/ii,/i 2 , . . . ,Hk}, respec- 
tively, As long as n > k and there are no mass degeneracies in the denominator, Z^ n,k ^ can 
be written as the sum of simple poles times their residues: 

^(WiWl-S^y, (32) 

j= l 1 ~r m i 

where 

^-( {m },{,}) S g^g,:g . (33, 

If there is a double pole, the residues are modified. We need consider only the case of 
one double pole; let it occur at q 2 = —mj. We then have 



d I Etife 2 + ft) 



dmj \ n?=i(? 2 + m ^ 



(34) 



Here the product over j includes £, i.e., 1 < £ < n. We now expand the quantity inside of 
the derivative as a sum of single poles and take the derivative of the resulting expression. 
The result is 



(m ■ fmWutt - ^MiM) + V ^ {{m} ' M) (35) 
l dp (m e ,{m}M)- {q2 + m 2 )2 +X, {q2 + m 2 ) > ( 35 ) 
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£4 

3=1 



n,k] 



with 

Dff ] ({m}M) = -^Rf M ({m}M) . (36) 
Note that Z)j™ takes on a simple form for j ^ I: 

Dff = Rf +1M ({m}'-M) (j + i) , (37) 

where {m}' is just the set {m} with repeated: {m}' = {mi, . . . ,m^,m£, . . . ,m n }. For 
j = DjJ becomes quite complicated, with n + k terms due to the differentiation. We 
emphasize that these formulae are needed solely for performing the momentum integrals 
explicitly in the partially quenched case. In full QCD, there are no double poles. 
We now collect some identities satisfied by the residues: 

1 , n = k + 1 ; 

, n > k + 2 . 

E"=i m ) - Ej=i A, n = k+l- 

E<'% 2 H-l, n = k + 2 ] 

, n > k+3 . 

1 , n = k ; 

, n > k + 1 . 

m t + EJ=i m ) ~ ELi »l , n = k ; 
J2{D^)-R^ = {-l, n = Jfe+l; (38) 

, n>k+2. 

These identities are easily obtained by expanding both sides of Eq. (jH2| or (jH5j) for large g 2 . 

When performing the explicit evaluation, the following integrals are needed: 

11 s / - i^<'" 2 > • < 39 > 

^ - / tSi^ = f tSz - ™^ - • (40) 



,[n,k] 
' j,' 

3=1 



(2vr) 4 g 2 + m 2 J (2tt) 4 16tt 
d 4 g 1 _ d 1 
(2tt) 4 (g 2 + m 2 ) 2 ~ ~dm^ 1 ~" 16~7r 

* - / (gi (?tW - - - ^ - , (42) 

where we have defined the chiral logarithm functions 

2 

£(m 2 ) = m 2 In — [infinite volume] , (43) 



fm 2 ) 



(Tfl 2 \ 
In — + 1 J [infinite volume] , (44) 
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with A the chiral scale. We use the arrow in Eqs. ()39j) through ()42|) and elsewhere to indicate 
that we are only keeping the chiral logarithm terms. If the system is in a finite (but large) 
spatial volume L 3 , the following modifications are required 14j : 

(Tfl 2 \ 
In — + 5i{mL) \ [finite spatial volume] , (45) 

(vr? \ 
In + 1 J + 5 3 (mL) [finite spatial volume] , (46) 



where 



(m 2 



. , T . 4 Ki(\r\mL) , Af7S 
5 3 (mL) = 2^K Q (\r\mL) , (48) 



with K and K\ the Bessel functions of imaginary argument. 
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V 



■d 
■d 



FIG. 1: The two-point mixing vertex coming from the U ' term, (a) corresponds to the chiral theory 
(we also have similar U — S and D — S mixing terms), (b) shows the corresponding quark level 
diagram. Here we only show the mixing among the taste- vectors, but there are similar vertices 
among the axial tastes, as well as the singlet tastes (with a 2 5' v — > 4ttt,q/3). 





(a) 



lb] 



FIG. 2: The SXPT diagrams contributing to the pion decay constant, coming from wave-function 
renormalization. The box represents the the axial current, (a) is the connected piece, where the 
propagator in the loop contains no two-point vertex insertions, (b) subsumes the graphs which 
have disconnected insertions within the loop. The cross represents one or more insertions of the 5' 
vertex, with 5' given in Eq. (JJJ). 



(a) 





FIG. 3: Same as Fig. |^1 but these contributions to the decay constant are from axial current 
corrections. 
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(e) 



FIG. 4: The quark level diagrams that contribute to the one-loop pion decay constant. The 
box represents an insertion of the axial current. The diagrams on the left correspond to the 
wavefunction renormalization while those on the right correspond to the current corrections. 




(b) 



FIG. 5: The quark level diagrams for 2 — > 2 meson scattering which contribute to f p +. The 
indices i and j represent arbitrary quark flavors. There are two additional diagrams (not shown), 
which are like those in (a) but have the roles of x and y interchanged. The box stands for the axial 
current. 
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